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Abstract 

Radiative heat transport between materials supporting surface-phonon polaritons is greatly en¬ 
hanced when the materials are placed at sub-wavelength separation as a result of the contribution 
of near-held surface modes. However, the enhancement is limited to small separations due to the 
evanescent decay of the surface waves. In this work, we propose and numerically demonstrate an 
active scheme to extract these modes to the far-held. Our approach exploits the monochromatic 
nature of near-held thermal radiation to drive a transition in a laser gain medium, which, when 
coupled with external optical pumping, allows the resonant surface mode to be emitted into the 
far-held. Our study demonstrates a new approach to manipulate thermal radiation that could hnd 
applications in thermal management. 
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Thermal radiation plays a role in many applications ranging from infrared detection and 
sensing applications for environmental and medical studies to energy harvesting with 
solar thermophotovaltics [SHS] and infrared emissions from Earth to space [Hj. Thermal 
radiation is also essential to thermal management applications as in microelectronics [7], 
space technology [8] and buildings [9]. 

In the far-held, the blackbody limit governs the maximum radiative hux between two 
bodies. Recently, a number of works have demonstrated that near-held radiative heat trans¬ 
fer is enhanced by many orders of magnitude compared to the far-held limit for closely 
spaced objects with either natural or engineered resonant surface modes [EHIE!. 

There have also been efforts to couple these near-held modes into the far-held with the use 
of grating structures na. antennas [TH], and a thermal extraction lens [Hj. However, all 
of the mentioned techniques to control radiation are passive and are still subject to the 
blackbody limit in the far-held. 

Unlike passive schemes, active schemes extract energy from a system through external 
work, allowing certain limits to be overcome without violating the second law of thermody- 
manics. For example, refrigerators move heat from cold to hot using electrical work input. 
In optics, external work in the form of laser light has been used to cool of gaseous mat¬ 
ter to sub-millikelvin temperatures [201 El] by removing kinetic energy from the atoms. In 
solid-state materials, optical irradiation can also cool materials by emission of upconverted 
fluorescence |22| due to removal of energy in the form of phonons. This concept has been ex¬ 
perimentally demonstrated to cool rare-earth doped glass jOOlEll to cryogenic temperatures 
and recently to cool semiconductors by 40 K from the ambient temperature izg. However, 
no active schemes are available to extract energy out of a system as thermal radiation. 

Here, we theoretically propose and numerically demonstrate an active thermal extraction 
scheme that extracts near-held thermal photons into the far-held. Our approach exploits the 
monochromatic nature of near-held thermal radiation to drive a transition in a laser gain 
medium, which, when coupled with external optical pumping, allows the resonant surface 
mode to be emitted into the far-held. Unlike the passive case, our active approach can lead 
to thermal radiative hux exceeding the blackbody limit in the far-held because work is being 
performed on the system. 

A schematic of the method is given in Fig.j^a). A laser gain medium containing emitters 
with discrete energy levels is placed in the near-held of a material that supports a resonant 
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surface wave. We model the emitters as a three-level system, as shown in Fig.[^b). An 
external pump laser is tuned to the 0-1 transition, exciting population into level 1 . If 
thermal radiation drives the transition from 1-2 and the 2-0 transition is radiative with high 
quantum efficiency, the electron transition will emit blue-shifted photons in the far-held, 
thereby extracting the trapped near-held thermal radiation. 

With a typical blackbody spectrum, the efficiency of such a scheme would be vanishingly 
small because of the low energy density and the broadband nature of thermal radiation [SE] • 
However, in the near-held, it has been demonstrated that the radiative energy density is 
nearly monochromatic and far exceeds that in the far-held by several orders of magnitude 
[ 27 j . Therefore, with near-held thermal radiation the 1-2 transition can be efficiently driven 
by matching the near-held energy resonance energy to the 1-2 transition energy. 

To study this system, we use rate equations to determine the steady-state populations in 
each energy level with external and near-held pumping: 


—— _ —tTL2(A^2 ~ Ni) — 7i2A^2 ~ 72oA^2 
^ = Wi 2 iN 2 - N,) - WoliN, - No) - 7ioA^i + 7i2iV2 
^ = Woiim - No) + 72oiV2 + 710A^i 

Nf = A^o + + A^2 


( 1 ) 

( 2 ) 

( 3 ) 

( 4 ) 


where IT12 is the absorption rate of the 1-2 transition as a result of the near-held energy 
density, IToi is the absorption rate of the 0-1 transition due to external pumping, W are 
population density of each level, Nt is the total population density for system and 7*^ is the 
overall (radiative and non-radiative) spontaneous decay rate of the i-j transition. Here, y’j 
stands for radiative rate of the i-j transition such that 7*^ = 7^- -|- We assume that all 
energy levels are non-degenerate so that Wij = Wji. Solving Eqs. to in steady state 
yields the equilibrium population densities for each level from which the power density can 
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Figure 1. (a) Schematic of the active thermal extraction scheme. An emitter with discrete energy 

levels is placed in the near-held region of a semi-inhnite planar substrate supporting a surface 
resonance. The external pumping couples with the near-held energy to be emitted as blue-shifted 
spontaneous emission in the far-held. Various other spontaneous decay channels into the near-held 
and the far-held are also taken into account, (b) Energy level diagram of the emitter for our 
proposed concept. The 0-1 transition absorbs external pump photons, and near-held photons drive 
the 1-2 transition. Spontaneous emission from the 2-0 transition emits near-held photon to the far- 
held. jij are spontaneous decay rates for the i-j transition with the superscript ”r” for radiative and 
”nr” for non-radiative rate. Wij are the absorption rate that depends on the energy density. The 
orange arrow indicates external optical pumping, the dashed arrows indicate various spontaneous 
decay channels with the blue arrows indicating the upconverted emitted photons carrying near-held 
energy into the far-held. 

be expressed as 

Poi = ^^ioW^oi(^o ~ ^i) 

_ (VFi 2(7 io + 72o) + 7io(7i2 + 72o)) 

^^12(710 + 720) + 710(720 + 712) + Wqi (3 VFi 2 + 2(720 + 712)) 

^20,net = ^(<^20 — ^10)720^2 

_ _ ^(^20 (^ 10 )NtWpi^20^12 _ /gx 

^^12(710 + 720) + 710(720 + 712) + Wqi (3 VFi 2 + 2(720 + 712)) 
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where Pqi is the external power density absorbed by the 0-1 transition and P 2 o,net is the net 
extracted power density into the far-held from the 2-0 transition. 

Using Eqs. |^and|^ the intrinsic efficiency of extraction can be expressed as the ratio of 
the amount of extracted energy radiated into the far-held by the 2-0 transition with respect 
to the external pump energy absorbed by the 0-1 transition 

_ P2O,net _ {^20 — ‘^10)720^12 _ /yN 

^ -Poi <^io(bUi2(72o + 7io) + 710(720 + 712)) 

In the ideal limit of a dominant radiative 2-0 transition 720 and strong near-held ab¬ 
sorption IU12, Eq. Intends towards (a;2o/<^io — l)(72o/72o) which depends intuitively on the 
ratio of the emitted net energy and absorbed photon energy and on the radiative rate of the 
2-0 transition for the photons that reach the far-held. When 7710 > 0, there is net energy 
extracted from the system if we assume no parasitic absorption of external pump energy. 
The intrinsic efficiency in Eq. depends only on internal parameters of the system and is 
independent of the absorption rate Wqi of the external pumping (0-1) transition. 

To estimate the efficiency of the scheme, we take properties based on rare-earth dopant 
embedded in gallium lanthanum sulhde (GLS) chalcogenide glass as the emitter system 
with typical values listed in Table [2HI121] • Such a system has a high quantum yield for 
transitions in the mid-infrared (MIR) region with a high melting point [HD] and has been 
used in making hber-based MIR lasers [2D]. Here, we assume the quantum efficiency of 
the equivalent 0-1 and 2-1 transitions to be unity and choose the wavelength-independent 
permittivity of the GLS chalcogenide glass to be 4.8 [31]. Then, we model the substrate 


Transition A(/ini) 

7°(s-') 

QE (%) 

0-1 

1.83 

1034 

100 

2-0 

1.22 

1370 

79 

1-2 

3.88 

36 

100 


Table 1. Parameters of a typical rare-earth emitter in GLS chalcogenide glass for modeling our 
proposed system. 7 )j(s“^) stands for the decay rate of the i-j transition for an isolated emitter and 
QE is the quantum efficiency of the transition. 

permittivity with the expression e{u)) = eoo{oj\ — LtP — i^oj) /(a;|. — uP — iju) where Coo = 5 . 3 , 
u!t = 388.4 X 10 ^^ = 559.3 x 10 ^^ s“^ and 7 = 0.9 x 10 ^^ s“^. We tailor the substrate 
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resonance to match the 1-2 transition with Re(es„6strate(<^)) = —^medium so as to enhance the 
energy density of the near-held thermal radiation with the emitter. 

To calculate the intrinsic extraction efficiency of this system using Eq.[^ we need to know 
near-held absorption rate W 12 , which in turn requires the near held energy of the substrate. 
To obtain this parameter, we use the formulation from Ref. [32] to calculate the near-held 
energy density I(aj) of the substrate at a temperature of 750 K where the blackbody spectrum 
peaks around 3.9 /rm, matching the 1-2 transition wavelength in Table|^ Fig. [^a) shows that 
the near-held energy density / (cj) increases by seven to eight orders of magnitude compared 
to the far-held energy density at very small distances from the substrate. Also, the spectrum 
peaks at a; = 485.5 x 10^^ s“^ when Re(esubstrate(<^)) = —^medium- Then, we can approximate 
the near-held absorption rate using of a recent formulation by Archambault et ah |33| which 
specihes the modihed spontaneous and stimulated transition rates by quantizing a surface 
plasmon for the planar interface. The isotropic stimulated rate in Eq. (29) of Ref. [3S| is 
also the absorption rate IT12 for the near-held in Eqs.[^to|^ We modify this expression to 
incorporate the energy per unit volume /(a;) = I{oi),k)dk in Fig.j^a) for the transition 
for different values of wave vector k to obtain 

k 




ij,near—field 


2hc(jQ 
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k)g{uj)dkd(jj 


( 8 ) 


(9) 


(u! — (Ac(j/2)^ 

where 7?- is the spontaneous decay rate for an isolated emitter and g{u) is the lineshape of 
the transition with a linewidth of Au). We follow the formulation laid out in Chance et al. 
to obtain the distance dependence of 7*^ of an isotropic emitter due to the modihcation of 
density of states by the surface in the near-held [HI] . 

To obtain the power densities using Eqs.[^and|^ we model the induced absorption rate 
lEoi due to far-held pumping using the well-known expression for the stimulated rate [35] 


Wij,external = ^ whcre 7^ is the radiative spontaneous decay rate that couples to 

external pumping from the far-held, A is the incident intensity of the external pumping 
held and n is the index of the chalcogenide medium. The linewidth for the 0-1 and 2-1 
transitions are assumed to be 2 x 10^^ s“^, comparable to typical laser gain mediums [35] . 

The intrinsic efficiency of thermal extraction versus distance from the emitter is shown 
in Fig.j^b). The maximum efficiency is small, around 0.33% and decreases to zero beyond 
a few hundred nanometers. The total extracted intensity is dehned as the integral of the 
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Figure 2. (a) Energy density at different distances d from the surface of the substrate. The top 

medium is GLS chalcogenide glass. The thermal radiation spectrum becomes nearly monochro¬ 
matic close to the substrate, consistent with other calculations [27]. (b) Extraction efficiency tjiq 
of external pumping from the 0-1 transition assuming properties in Table [I] The low efficiency in 
the blue line is a result of the enhanced spontaneous rate for 1-2 transition in Fig.[^ (c) Integrated 
power extracted for emitters uniformly distributed from 10 nm onward from surface. The density 
of emitters is assumed to be 10^^cm“^. The saturation behavior approaches the green dashed 
’’saturation” line due to the finite number of emitters in the system saturating the population 
difference at high input powers. 


power emitted by the 2-0 transition over all distances, P 20 ,netd^ = idl 2 o^net/dz)dz = 
^^20,net- We integrate from 10 nm onward until the intrinsic efficiency decreases to almost 
zero. Figure c) shows the extracted power per unit area as a function of inpnt power /„. 
The extracted power increases linearly with the inpnt power for low power inpnts before 
satnrating at higher powers, but the overall power extracted is orders of magnitude lower 
than the input power. The saturation is due to a high pump power which causes the 
population difference between the levels to tend to zero, leading to a maximum in absorption 
and thns the satnration behavior. A limiting case of Eq.|^ can be fonnd for large VFoi as 
h{u2o — <^io)bFi272o-^t/(3bFi2 + 2(720 + 712))- Integrating this limit over distance agrees with 
the satnration as plotted in Fig.j^c). 

Figureshows that active thermal extraction is possible, but both the intrinsic efficiency 
and the total power extracted are very small for the chosen parameters. However, according 
to the limit of Eq.[^ the maximnm efficiency shonld be aronnd 35%, much higher than in 
the example. To nnderstand the reason for this difference, we examine Eq.j^in more detail. 
The maximum efficiency occurs when 720 and IF12 are large. We calculate the transition 
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rates versus distance from the substrate in Fig.|^a), and observe that the transition rates 
for 0-1 and 2-0 transitions are not affected by the presence of a surface as they are off- 
resonant. However, the decay rate for the 1-2 transition 712 is strongly enhanced as the 
emitter approaches the surface [SH |36l |37j. As a result, the near-field absorption rate is 
smaller by about two orders of magnitude compared to the decay rate even though both are 
enhanced by orders of magnitude due to the increase in the optical density of states in the 
near-field. Physically, this calculation indicates that as electrons are excited from energy 
level 1 to 2, they immediately decay back to level 1 at the rate 712. 

The reason for this cycling is that the thermal near-field energy density is not sufficient 
to allow near-field absorption to dominate over near-field spontaneous decay. Archambault 
et ah [32] also highlight the need for some minimum energy density for stimulated emission 
to dominate spontaneous decay. Unlike the case for stimulated emission of surface plasmons 
with external pumping such as in Ref. [3H| where the external laser field intensities can be 
tuned, here the thermal energy density is restricted to that for a blackbody. Thus, the 
spontaneous decay rate will always dominate over near-field absorption for realistic values 
of near-field energy density. On the other hand, Fig.|^a) also shows that while a resonantly 
enhanced 712 offsets the enhanced absorption VF12, the extraction efficiency 7710 still requires 
a large value of lUi2- Beyond a emitter-substrate distance of about 100 nm, the extraction 
efficiency in Fig.j^a) drops significantly as a result of the low near-field energy density, 
although the ratio IU12/712 remains in the same order of magnitude up to a distance of 1 
/rm. 

Therefore, to break the cycling between levels 1 and 2, it is essential that the strongly 
radiative decay rate from 2-0 (720) is comparable to the decay rate 712 in the near-field. 
Figure l^b) shows that the efficiency is boosted to almost the ideal limit at short distances 
if 720 is increased substantially. In Eq.[^ if we increase 720 to be more comparable to 712 in 
the near-field, then the ratio of 720/720 begins to dominate in the expression, increasing the 
extraction efficiency towards the ideal limit discussed earlier. 

The factors discussed above affect the intrinsic efficiency, but the total extracted power 
also depends on the input power Hoi and the emitter density Nf. We now examine the 
role of these parameters. Firstly, the absorption of the pump power Hqi depends on the 
linewidth of the 0-1 transition, and decreasing the linewidth increases Hqi in Eq.j^due to the 
increased concentration of input power in a given bandwidth for each emitter. Figure |^c) 




Figure 3. (a) Normalized spontaneous decay rates versus distance for three different transitions. 

The 2-1 transition is on resonance with the substrate dispersion and is enhanced greatly whereas the 
0-1 and 2-0 transitions are not significantly affected by the presence of the substrate, (b) Intrinsic 
extraction efficiency 7710 versus the scaling of the spontaneous rate 720 at d = 20 nm. The blue line 
shows real behavior according to Eq.[^ Increasing 720 greatly enhances the efficiency so much so 
that around 10^720 the efficiency approaches the ideal limit of the system, (c) Integrated extracted 
power versus the linewidth of the 0-1 transition for a input power of 0.01 W/cm^ in the linear regime 
with the same parameters as Fig.j^b). Decreasing the linewidth increases the absorption efficiency 
and therefore the output power, (d) Integrated extracted power for different emitter densities in 
the saturation regime with an input power of 10^ W/cm^. The full description using Eq.|^at high 
input power agrees with the saturation limit shown in the dashed green line. Increasing the emitter 
density increases the saturation extracted power by the same order of magnitude. 

shows that decreasing the linewidth of the 0-1 transition increases the integrated extracted 
power as predicted. Secondly, the total dopant density Nt also affects the extracted power. 
As discussed earlier, the saturation limit at higher incident powers is proportional to the 
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Figure 4. (a) Intrinsic extraction efficiency versus emitter-substrate distance for an optimized 

system. The extraction efficiency follows the ideal limit for small distances before decreasing due 
to a decreasing Wi 2 > The extraction efficiency is much improved compared to Fig.j^b). (b) 
Integrated power extracted of the optimized system with emitters uniformly distributed from 10 
nm onward from the substrate surface. A reduced 0-1 transition linewidth of Awio = 2 x 10® s“^ 
and a higher emitter density Nf = 10^^ cm“^ lead to much higher saturation limit shown in the 
dashed line. 

dopant density, and therefore the dopant density must increase to increase the saturation 
limit. Figure|^d) shows the integrated extracted power at a high input power as a function of 
emitter density. The extracted power scales accordingly with the dopant density and agrees 
well with the prediction from the expression for the saturation limit. Thus, increasing the 
dopant density shifts the saturation limit in Fig.j^c) up proportionally. 

Using this understanding, we now recalculate the efficiency and extracted power for an 
optimized gain medium with the spontaneous rate for the 2-0 transition increased to 1.37 x 
10® s“^, Ac(Jio = 2 X 10® s ^ and Nt = 10^^ cm Figure a) shows that the intrinsic 
extraction efficiency is much higher than in Fig.j^b) and almost near the ideal limit for 
small emitter-substrate distances. The decrease of efficiency at larger emitter-substrate 
distances is due to a decrease in near-field coupling. Figure |^b) shows a much-increased 
integrated extracted power at each given input power compared to Fig.j^c). The saturation 
limit derived earlier also agrees with the full calculation at higher input powers. This limit 
is independent of the input power and is proportional to the dopant density. 

This calculation shows that the active thermal extraction scheme has potential to ef- 
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ficiently extract a significant amount of near-field thermal radiative energy. The key to 
realizing this potential is to identify an appropriate gain medium with surface resonance 
and a emitter with matching transitions in the mid-infrared wavelength range where pho¬ 
tons are thermally populated at typical temperatures. Decreasing the linewidth (such as 
using a cavity) or increasing the 0-1 spontaneous decay rate helps to decrease the input 
power necessary for extracting a given power, but a high dopant density is still required 
to increase the saturation limit. Cerium doped crystals can potentially be a candidate as 
they have a 4/°5d^ 4/^5d° transition with a short lifetime of around 40 ns [HH], ideal for 

the 2-0 transition proposed here, as well as a mid-infrared transition of 4.5 jim (HH] for the 
near-field absorption. Radium is a very good candidate with parameters such as transition 
wavelengths and decays rates comparable to our optimized model. It has been recently 
demonstrated to interact with black body radiation in a laser trapping experiment HJ]. 

Our work shares some similarities with laser cooling of solids [2SH2SI SI] and active 
schemes in plasmonics photonic crystals [l3|, and metamaterials mi SS] but differs 
in a number of important ways. First, laser cooling directly extracts phonons, while our 
scheme extracts surface phonon polaritons. Therefore, our scheme has potential to be much 
more efficient than laser cooling because of the significantly higher energy of surface phonon 
polaritons than phonons. Also, laser cooling requires the medium to be cooled to possess 
very specific energy levels, whereas our scheme only requires that the medium possess a 
surface resonance. Second, active schemes in plasmonics have been used to realize spasers 
and to compensate loss. One main distinction between these schemes and our approach is 
that here one transition is pumped by a near-field, incoherent thermal radiative source rather 
than a coherent pump. As a result, our approach does not lead to any form of stimulated 
emission or coherent single mode emission. Additionally, the thermal energy density is not 
sufficient to cause the imaginary part of permittivity of the gain medium to become positive; 
our medium is actually absorptive under all conditions. 

In conclusion, we have numerically demonstrated an active thermal extraction scheme 
that allows bound surface waves to be converted from evanescent to propagating waves. 
Our work exploits the monochromatic nature of near-field radiation to drive a transition 
in a gain medium simultaneously with an external pump. This demonstrates shows the 
high potential for manipulating thermal radiation using active processes rather than the 
traditional passive approaches. 
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